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ABSTRACT 


We  give  rational  approximation  schemes  interpolating  the  values 


of  a given  function  along  two  edges  of  a triangle  and  the  normal 


derivative  along  the  third  edge.  In  addition,  we  give  error  bounds  for 


our  schemes.  For  the  uniform  norm,  these  bounds  are  best  possible. 
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BLENDING  INTERPOLATION  SCHEMES  ON  TRIANGLES  WTTH  ERROR  BOUNDS 


K.  Bohmer  and  Gh.  Coman 


0.  Introduction: 

A feu'  ye-jrs  ago  Barnhill-BirKhoff-Gordon  [1]  constructed  rational 
functions,  interpolating  the  values  of  a given  function  on  the 
boundary  of  a certain  normalized  triangle  T.  In  the  same  paper  inter- 
polation schemes  are  given,  realizing  the  values  of  a given  function 
and  its  normal  derivatives  on  the  boundary  of  T.  These  schemes  are 
affinely  invariant  iff  the  normal  derivatives  are  transformed  into 
the  corresponding  directional  derivatives. 

/ 

Here  we  give  rational  approximation  schemes  interpolating  the  values  of 
a given  function  along  two  edges  of  a triangle  T and  the  normal 
derivative  along  the  third  edge.  In  addition  we  study  the  remainders 
of  our  schemes  and  give  better  error  bounds  than  in  [6].  For  th.e  L^- 
norm  the  error  bounds  are  best  possible. 

2 

As  triangle  T we  choose  the  standard  triangle  Tj^={(x,y)6  IR  |x  > o, 
y > o,  X + y < h}  with  vertices  (o,o),  (o,h)  and  (h,o).  Our  schemes 
are  affinely  invariant  in  the  same  senee  as  discussed  above. 


1.  First  Interpolationscheme ; ■ 

In  contrary  to  the  following  problems  we  treat  Problem  1 fairly 
extensively . 

Problem  1 : For  f : IR  uith  exieting  on  x = o find  a 

blending  function  Gjf  such  that  ^ 


(1) 


1 


Vly=o  = ^>y=o  ® * x=o 

V'xty=h  = "'x.y^h^""  " " • 
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AX 


As  in  many  similar  cases  we  use  the  Boolean  sum  of  two  operators 
namely 

a : J J 


a® (X  GLis  , 


X : j * j- 


Here  we  take  CL = Pj , = P2  with 

^l^'x+y-h  ■ ^U+y  = h ' ^2^lx+y  = h ’ 

•x  = o 'x=o  ’ 

*’2^ly=o  ■ ^ly=o 

Immediately  one  finds  (compare  the  technique  used  in  [1]) 
(Pjf)(x,y)  = f(h-y,y)  + (x+y-h)f (o ,y) 
(PjfKx.y)  = f(x,o)  + f(x,h-x) 

and  therefore 

C-jf  = Pjf  © ^2^- 

We  find  for  Gj^f  with 

.1 

Uj  :=  Pjf,  Uj  :=  Pjf 
the  formula 

(PjPjfXx.y)  = (PjUjXx.y) 

= f(h-y,y)  + (x+y-h)  j ^(f (o ,h)-f (o ,0) ) 


^ f^^‘°\o,o)  + ^ (f^^‘°^(o,h)  - f^°’^\o,h))l 


and  so  finally 


(GjfXx.y)  = ll^^f(x,o)  ♦ j^f(x,h-x)  ♦ (x+y-h)f *°\o ,y) 

+ (h-x-y)|  ^(f(o.h)  - f(o,o))  4 »i^f<l*o>(o,o) 

♦ J(f<l.o>(o,h)-f^°»^\o,h))] 


J 


Now  it  is  straightforward  to  check  that  safisfies  (1)  and  it  is 
equally  simple  to  prove  the  rest  of  the  following. 


Theorem  1 : l^et  f : Tj^  •»  IR  he  given, 

(a)  If  f^^’°^(o,y)  for  y € lo,h]  and  f^°*^^(o,h)  exist,  then  G.f 

in  (2)  patisfies  (1).  "orcover,  if  f(*,o),  f(*,h-*)  and  f^  ’°\o,‘) 
e C^[o,h]  then  G^f  e C(T^)  for  j = o and  G^f  6 (T^\  { (h  ,o) ) ) 

for  j > o.  •« 

(b)  G reproduces  bivariate  polynomials  of  total  degree  at  most  2. 

• 1 

(c)  f/  f (•  ,o)  , f (•  ,h-’  1 f ’°No,-  ) 6 C^[o,h]  , then  is  a 

solution  of 

2 

4,  r o tn  T \ {(h,o))  . 

X y 

(d)  If  f(',o),  f(o,‘),  f^^’°\o,*)  6 L^[o,h)  then  G^  satisfies  the 
following  kind  of  a maximum  principle 

U(Tj^)  = II  'I  Uolo.h]  * 11^^°’*  ^ II  U,[o,h] 


♦ hll  f^^’°^(o,0  + ^f^^’°^(o,o)ll 


U,[o,h] 


^1  f (o,h)-f (o.o)  + h(f^^’°^(o,h)  - f^‘^>°^(o,h))l  • 


3.  Remainder  for  the  First  Interpolation  Scheme: 

To  estimate  the  remainder  R^f  in 

(3)  Rjf  :=  f - Gjf  for  G^f  in  (2) 

we  use  a ''Sard-Kernel-Theorem"  [ 16]  on  triangles  due  to  Barnhill- 
Mansfield  [2]  : Let  F be  a linear  functional  of  the  form 

(4)  F(g)  = I a g^^^^x,y)du  . . (x.y)  + 

i<  p T 

j<  q 

♦ I f (x)  + I I g^^*3^a.y)cV:.(y), 

i + j<m  I.  ■'  i*j<m  I,  ■' 


where  the  functions  are  of  bounded  variation,  T is  a given 

triangle,  (a,b)  6 T is  such  that  the  rectangle  with  opposite  vertices 
<a,b)  and  (x,y)  and  with  edges  parallel  to  the  coordinate  axes  is 
contained  in  T for  all  (x,y)  6 and  1^ , I2  are  the  intersections  of 
the  line  y = b respectively  x = a with  T.  If  F satisfies 

F(p)  = 0,  p € = (polynomials  of  total  degree  ^ m-1)  > 

then  we  have 

C5)  F(g)  = I J l^"^’^(x,y;t)g^'""^»^^(t,b)dt 
j<q 

I * 

, ♦ I / K^’"*"^(x,y;T)g^^»"'"^\a,T)dT 

i<p  Ij 

«. 

' ♦ |/  KP‘‘l(x,yit,T)g^P‘‘l^t,T)dt  dT,  g e 


where 


fi 


We,  too,  will  use  the  class  B^’*"  = B!'’'"<a;b)  of  functions,  r > 1, 

p.q  p*q 

whith  slightly  modified  properties  (7)  namely 


g 0 < j < q 

g(i,m-i)^a,.)  e L^dj),  0 < i < p 


€ L_(T). 


Now  we  are  able  to  prove  the  following 
Theorem  2 : f/  f 6 Bj’2  then 


(8)  llRjfll  , , , 5 5;i|  ^ = 

«*  n ^ «o  n ^ 


L [o,h] , 


2 3 

where  the  constanta  are  best  possible.  If  f t then 

(9)  llRifll  L (T  ) = L (T  ) ■" 

^ 6*/3o  2'  h'^ 


Ljlo  ,h] 


Remark  1 : Best  possible,  but  unusual,  estimations  for  the  Lj-norm  are 
valid  if  f 6 B^’^Cojo): 


’I 


2,1 


’L„[o,h]- 


*■1*2  *1  \ 

This  is  proved  by  (13),  (17)  via  JJ(g(x,y))  dxdy  with  g * from  (15) 
and  similar  arguments  for  the  other  indices. 


-5- 


Proof : Since,  by  (2)  and  (3),  R^f  is  of  the  form  (4)  with  m = 3,  p 
q = 1,  a=b=o  we  find,  using  (5), 

f (Rjf)(x,y)  = / K^jJ°(x,y;t)f^^»°^(t,o)dt 


(10)  -I  ♦ / K°»^(x,y;t)f^°*^^o,T)dT  + / (x  ,y  ;t  )f  ^ (o  ,t  )dt 


Now  (6)  implies 


K^h’°<x,y;t)  = R 


[^•] 


(x-t)»  _ h-x-y  (x-t)|  _ y (x-t)^  _ 

2 h-x  2 h-x  2 


^*3(x,y;t)  = R^j^<J^+j  = IVJLl*  - { 


y (h-x-x); 
h-x  2 


(h-x-y)y (h^-T^) 


(11)  ^ KV^(x,y;x)  = R^(x(y-x)^)  = x(y-x)^  - (h-x-x)^ 


•»•  (x+y-h)  (y-x)^  + (h-x-y)^  (h-x)^  | 


= (y-T)^(h-y)  - (h-x-x)^  - J (h-x-y)(h-x)  , 


I KV^(x,y;t,x)  = Rj((x-t)^  (y-T)°)  = (x-t^[  (y-x )°-jj^(h-x-x )°] 

Since  the  proof  is  a little  bit  lengthy  we  want  to  scetch  what  we  are 
going  to  do:  (11)  implies  that  there  are  three  nontrivial  summands 
in  (10),  which  we  discuss  separately.  We  first  show  how  to  find  sharp 
error  bounds  and  then  treat  the  kernels  , K°j\^  and 


I 


Sharp  error  bounds : To  find  sharp  er'ror  bounds  we  have  to  look  at 

(RjfHx.y)  = (R°>^f)(x,y)  + (R^’^f)(x,y)  + (R^’^f)(x.y) 

with  (R^’^f)(x,y)  defined  corresponding  to  in  (10)  and  (11).  Let 

for  example 


T^:=T^^«l.-T2;|^^y):={(t,T)  e T J ^ ( x ,y  ;t , t ) j 


2.1, 


ol 


T_:=T^>^;=T^’J^^yj  : = {(t,T)  6 (x  ,y  ;t  ,t ) < o}  . 


Then 


I (R^’"^f )(x,y)|  <//  + //  |K^d^(x,y;t,T>f^^’^^(t,T)|  dt  dr 

- .j,  ,p 


- **  ) * -^^  I K^j;^(x,y;t,t)  | dt  dt 

^ h T^ 


or 


(12)  Hr  ’ ’^"b„{T^)  ' II  I ^^h  .T )1  dt  dt 


since  f 6 the  error  bound  (12)  is  sharp.  To  estimate  ||  ^ f || 

1 > c ^2  ^ . 

we  look  at 

2 2 
((R^’^f)(x,y))  = (//  K^’^(x,y;t,t)f^^’^\t,T)dt  dt) 


.2.1, 


< //(K^’^(x,y;t,t))^dtdt  • llf^^’^^ll  ^ , . 

■ Tj^  ^2^h^ 

or 


Here  the  first  inequality  is  sharp  for  every  fixed  (x,y)  € Tj^ , the 
second  for  all  (x,y)  e T^.  So  we  finally  have 


<13) 


o ^ 2 1/2  -V 

<(//(//  KY(x.y.t.T)dtdx)  dxdy)  • IK  ’ IIl  (T  ) 
- ~ ” •>  h 


Th 


or 


- JJ  (/J  (K^^^(x,y,t,T))  dtdi)dxdy 

'*  «wt  m *• 


1/2 


Th 


[f(2.1)|, 


In  the  first  line  the  constant  is  best  possible,  whereas  it  is  not 
in  the  second  line. 

Similar  arguments  hold  for  R°’^f  and 


Now  for  0 

< X , o ^ 

0 

HA 

X 

1 

x: 

0 

liA 

X 1 

1 

|.o , 3 
h 

'H’h’H’  ■ 

h ^ K°j^(x,y;x) 

and 

_ V®  » ^ 1 

■ ^ h *h=l 

^1,2 
^ h 

fX  y .T  . . 

^K’n’n^  - 

h"^  K^j’^(x,y,x) 

and 

■ ^ h 'h=l 

^2,1 

h 

/X  y .t  T , 

= h"^  K^j’^(x,y;t,x 

) 

and 

^2,1  . „2,1 
^ 1 ■ h 

1 IK^^’^^IIl  (T  ) ' I I dtdT 

h Ho  h T, 


So,  for  example,  by  (12) 
and  with 

;=  {(x,y)  e IR^  I o < x,  o < y,  o < l-x-y>  = , 


X*  • “ • V*  ‘=5^  t ' • = « T*  :=  r 

-’R’y  • h’^  • H‘  K 


we  have 

fl  1 K^^^(x,y,t,T)  1 dtdt  = 

’’h 

= ^/  I K^’^(x’ ,y' ;t' ,t’)  ) hdt’hdt 


= //lKV^(x’,y';t’,T*)|dt'dT'  . 

rp  1 


-8- 


In  a similar  way  we  have 


’’h  ^ 


and 


//^//<X°h^(x,y;T))  dDdxdy  = 


■^h 


h“  //(//(K°’^(x',y';T'))  hdT ’ )hdx 'hdy ' 

fT>  fT*  •• 


= h’  ;/(//(K°’^(x' ,y' ;t'))  dT')dx'dy’  . 

m rn  ^ 


T,  T, 


So,  to  get  the  constants  in  (8)  and  (9),  it  is  enough,  to  confine 


So 

X y 1 *"  X 

niK^’^(x,yit,T)ldtdT  = / < / i^^(x-t)dT  + / •r^(x-t)dT  )dt 

Tj  t = 0 T=0  ^ Try 


LIXJIZ  V (x-t)^l’' 

1-x  ^ ‘t=o  1 - X 


x^y (1-x-y ) 


and  therefore 


i <T,.  <x.yi  e T.  . 


To  find  llg  * 11^  (X  ) 


CIS)  g^’^(x,y):=  X yCl-x^v)  ^ ^ 6 , = o iff  (x;v)e  3 , 


we  derive: 


2 1 2 

i&_  = ^ — (l-x-2y)  = o iff  X = o or  1 - X : 2y 
3y  1-x 


— 5-  [2(l-x)^  + y(x-2)]  = o . 
(1-x) 


2 1 2 1 

So  (x,y)  = — (x,y)  = o iff  x = o or  8y^-2y^-y  = o, 


1-x  = 2y  , 


12  • 2 1 
s.t.  y=-g,  x=-jis  the  maximum  point  for  g ’ , so 

= ^and 

«o  J 


(16)  l|f”'‘)lL.(Tj) 


To  find  estimations  for  llR^’*fl|i  (t  ■)  '~a  proceed  as  follows  (see  (lU)) 

2^  1 


[ (R^’^f)(x,y)]  ^ ^ K^;^(x,y;t,T)f^2’^^(t,T)dtdT) 

t=0  T=0  ^ 


X 1-x  o i ^ ^ i A > (j  ^ ^ 

< / / [ Cx  ,y  ;t  ,T )]  dtdt  / / [f'  * (t,T)5dtdt 


X 1-x 


t=0  T=0 


t=0  T=0 


1 / ( f ^(x-t)^dt  ♦ (j^)^(x-t)^dT)  dt  llf^^’^^i 


t=0  T=0 


yC)-x-y)  Ilf  ( 2,1)  „ 2 


3(l-x) 


SCTj) 


So  finally  HR  ’ flii  (t  ) = 


^1 


LjtT^)  ' 18. 6o' 


or,  with  g 


2.1  from  (15)  and  //  (g^ (x ,y ) ) ^dxdy  = 1/Uoo.21 

T. 


llp(2,l) 


>4(T^) 


(17)  1 


:(2 ,1) 


2o/TT 


ad  ^ : From  (11)  we  have  for  h = 1 

(18)  K°^’^(x,y,t)  = 2^ - 2(l-x")  ^ 2^ 


We  first  prove  that 

(19)  K‘^j^(x,y;T)  > o for  (x,y)  6 and  t 6 [o,l]  . 


In  o < T < y we  have 


K°i>^x,y;t)  = - jfl-r^ry  (1-x-T)  + 


r(l-X-v) (1-f 
2 


<°i’^(xiy  ;t  ) 


is  linear  in  t 


3 K°i’  ^ ( X , y ; T ) 


3K°i»^x,y-,T)  , 1 


= y(l-x-y)(y!-  - y)  > o in  Tj 

T=y  1"X  - •* 


and  = o iff  y-o  or  1-x-y 


K°’^(x,y  ;t  ), 


r 

This  implies  K°’^(x,y;t)  > o in  o < t < y. 


For  o<y<T<l-xwe  have  (see  (18>) 

‘ K°^^(x,yjT)  = - ■nfr3n‘  - (l-x-T+T-y)(l-T^)(l-x)] 

S + „ ^y_^-y  [(1-X-t)(1-T  + Xt)t  + (T-y)(l-T^)(l-x)l  > O . 

For  o<y<l-x<T<l  K°’^  reduces  to 

K°.3u,y..,)  = ysir^zilllzjll  > o.  • • 

so  (19)  is  proved. 

Since  K°’^  is  a semidefinite  kernel  we  have  with  n 6 (o,l) 


(R°’ 


^f)(x,y)  = / K°^^(x,y,T)f^°’^\o,T)dT 


= f^°*^^(o,n)  / K°j*^(x,y;T  )dT 

0 

1 

and,  analogously  to  (12)  and  (13), 

(X  |)f^°’^\o,*  )llj^  (o,hj'l  ^ ) 

m n «o  o ” n 

and 

Since  we  are  interested  only  in  the  constants,  it  is  enough  to  discu 
h=l  in  Corresponding  to 

K°i’^(x,y;t)  = (4>j  + ♦j  + <fr3)(x,y;T)  with 

2 

♦j(x,y;T)  :=  = o for  y < ^ 1 1 

= o for  i - x < t < i 


. v(l-x-v)(l-T^) 

♦ o(x,y;t):=  2 2 


we  find 


\ ^ Y 

G°’^(x,y)  :=  / (K°J^(x,y,T))  dr  = / * 2^i«3>dT 


o 
1-x 


♦ / (412  + 2({)2({>3)dT  ♦ / (Jij  dx 


- y - _y 


1-x 
6 


■ |^2o  ^ 12^  ^ ■ 6 6o(l-xy 


+ (1  - 

6 lo 


y (1-x)  y (1-x-y) (1-x)  (1-x)  . 

2o  ' 6 lo 


+ 2y  (1-x-y ) 
lb 


and  finally 

ff„o,3,  vjj  397 

//  G (x,y)dxdy  = 

1 

(20)  4 


/ 


397 


317S200 


/397/2 

12bo 


Computing 


G°’^(x,y)  :=  / K°’^(x,y;T)dT  = f (4>^  + <J> 


-2  ^ 


1 

(if,  + if,)dT  + / If  dx 

^ ^ 1-x  ^ 


we  find 


= y(y^  - (1-x)^  + 2(l-x-y))  = ^[(y-1)^  " ] > o 


S°’^x,y) 


2 2 

since  1-x-y  >o  and  therefore  (y-l)  > x . Because  of 


(x,y)  - Z (+  2(1-x)-2)  = - ^ = o iff  x = o or  y = 
3x  b J 


the  extrema  of  in  are  realized  on  3?^: 

- 9(0,3 , , _ y(y-l)^  ^ 2 

o < G ’ (x  = o,y)  = i gj  , 

^°’^(x,y=o)  H o , 

&°*^(x=l-y ,y)  = o . 


So 


(21)  llS°’^(-,v 


) * 81 

«o  1 


ad  With  h = 1 (11)  implies 

(22)  K\’^(x,y,T)  = (y-T)^  (1-y)  - ^ (1-x-t)^  - (1-x-y)  y(l-T ) 


Since  is  a continuous  piecewise  linear  function  (for  x ^ 1), 

one  finds 


(23)  K\’^(x,yit)  < o 


1 2 

by  simply  checking  K ’ (x,y  ; t ) for  x = o,  y,  1-x,  1. 


Again  we  have  to  compute 

19  1 1 , 2 

^^’^'(x,y)  :=  f 
o 

1 


G^*^(x,y)  :=  / (x ,y ;t ) ) dt  and 

o 

6^’^(x,y)  :=  - J K^^^(x,y;x)dT  (see  (23)) 


and  find 


G^’^(x,y) 


= ^ [<•- 


y)y^ (-2  + 3y-xy-y^ 


. ^ /19/14  ^ 

(24)  //  G^’'^(x,y)dxdy  = ( 2o  ) , 

T. 


♦ xy^(l-x)(2 -2y-x^-xy>  ♦ (l-x-y)^y^j 


o . 


* J(x^-y^*2y-l)  = ^(x*-(l-y)^)  | 

Because  of  aS^’^/3x  s o iff  x = o or  y * o we  find  j 

by  discussing 

o < s^*^(x=o,y)  = ♦ = * It 

^^’^(x,y=o)  i o 
*^(x=l-y  ,y)  = o 


or 

(25)  11^^’^IIl  (t  ) = It  • 

«o  1 


By  fitting  together  formulas  (16),  (21)  and  (25)  we  have  (8),  by  (111, 
(20),  (24)  we  have  (9)  and  the  constants  are  best  possible  in  the 
sense  indicated  in  Theorem  2 and  in  the  following  Remark.  □ 


4.  Second  Interoolation  Scheme: 


By  interchanging  x and  y we  come  from  Problem  1 to 

Problem  2 : For  f : IR  uith  existing  on  y-o  find  a 

blending  function  G2f  such  that 


vix=o  = ^ix=o  y ® 


(26)  ^ ®2^lx+y=h  ' ^lx+y=h 


(G,f)^°»^^l 


y-o 


'y=o 


The  corresponding  results  to  Theorem  1 and  2 are  valid  , too, 
and  we  have  for  Gj  : 


(Gjf)(x,y)  = f(o,y)  + ^ f(h-y,y)  + (x+y-h)f^°*^^x,o) 


h-y 


(27)  < 


r w / X (o,i)  (i,o) 

♦ (h-x-y)  U^(f(h,o)  - f(o,o))  + f^°’^'(o,o)*^(f (h,o)  + f (h.o)) 


-15- 


5.  Third  Interpolation  Scheme  and  Remainder; 

How  we  discuss 

Problem  3:  For  f ; ♦ IR  uith  oxiating  on  h-x-y  = o 

n y 


find  a blending  function  G^f  with 


S^U--o  = ^'x  = o ^ ® 


(28) 


x=o 


Ggfly,^  = fly,o  for  y e to.h] 


\ 


\ x+y=h 

V 

\ 

V 


y=o 


(G, 


^ ^'h-x-y=o  - 'h-x-y=o- 


h-x-y: 


By  distinguishing  the  two  different  cases  x ^ y and  x < y we  define 
Gjf  and  Gjf  by 

Olf|y=o  = fly=o' 


Ql^lyrO  “ ^ly=0* 

* *^2^  'lh-x-y=o  ^ Vx-y=o’ 

and  we  find 

^ _v,  ^ X ■ ../y^(l.o)  . ^(o.l)  X ,h+x-y  h-x+y^ 

(Qj^fXx.y)  = f(x-y,o)  + y(f  ’ + f ’ )( — , 


(29)  < 


X X J.  pto>l)x  /h  + X-y  h-X  + Vx 

(Q_f)(x,y)  = f(o,y-x)  + x(f  + f ’ ) ( — 5—^, — 5^) 


and  finally 


(30)  (G3f)(x,y)  :■ 


Qjf(x,y)  for  y < x , 
Q2f(x,y)  for  y > x . 


-If- 


A straightforward  computation  yields 

Theorem  3:  Let  f : T.  ♦ ]R  be  given  and  ||f  ( • ,h- • ) ||,  r . , =supess|  f (x  ,h-> 
7,  (Ax  x6[o,h] 

(a)  Jf  (f'  ’ + f ’ exists  along  h-x-y=o, then  G^f  is  a solution 

of  Problem  3.  Moreover,  if  f(o,*)  and 

(f(l.o)  ^ ,h-* ) 6 C^[o,h],  then  G^f  6 C^(T^). 

(b)  reproduces  linear  functions. 

(c)  Jf  f(*,o),  f(o,*),  ♦ f ,h-*  ) e L^[.o,h]tbcrj 


+ h|((f^^*°^  + f ,h-- )|| 


L lo,h] 


To  discuss  the  error  we,  again,  introduce 


(31)  Rgf  :=  f - Ggf  . 

Since  we  have  to  transform  the  original  triangle  Tj^,  we  have  to 

transform  the  conditions  defining  B^’^(a,b),  too.  Let 

P >4 

f:T^-IRlo<i,j: 

f 

C(Tj^)  for  (i,j)  < (p,q)  % 
AC(Ty^  {x  = yl)  for  j<q, 

i j 

(32)  B^:::  ==  ^ ^ ^ tx=y))  for  j<q,  i=m-j  , 


AC(T.  {x+y=h})  for  i<p,  j=m-i-l, 


^r^^h  i<P*  i=m-i. 

= (p,q) 

Theorem  U:  Jf  f 6 Bj  * then,  with  [jf  ( • ,h- • ) |L  . , ,,=supess  |f(x,h-x)|, 

2 ^ x6[o,h] 

I llR3f  lll^(T  ) % F""^Tx  ^*'u»[o,h] 


♦ ^iy  ' k^^"L,(T.  ) . 


3 


-17- 


1/  f e then 


(34)  i 


a 3 ^ Ml 

^ » (1-  * I-)  f(-,h-->ll 


'Ljlo.h] 


• All<k*l7>'<ly-S>^"4<v'' 


48/5 


The  constants  given  in  (33)  arid  (34)  are  best  possible  and  a sinilar 
remark  to  Remark  1 is  possible. 


Proof:  To  be  able  to  apply  (4)  and  (5)  we  transform  the  variables: 


x+y  ..  y-x  ^ - u^v  . u+v 
(35)  » - 7f  • ^ - /j  ■ /2 


Transforming  Q,f  via  (35)  wo  find  with  (29).  (30)  for  y > X or  v , o 


(“(jjfXu.v)  = g(v,v)  ♦ (u-v)g*^’°h^,v) 


with 


\ 11“  V U + V \ 

g(u,v)  :=  f(x  = — , y = 


(36)  < 


((J,f)(u=v,v)  = g(v.v)  = f(o,/2v)  for  v 6 [ o 


s 1_  + f^°’^^)(x=h-y,y)  . 

/7 


Applying  3*^*^  (5)  to 


(1i3?l2f)(u,v)  :=  (g  - 


= g(u,v)  - g(v,v)  - (u-v) 


with  m = 2 = p.  q = 1.  a = and  b = o we  have 


(?(3Q2f)(u,v)  * / Kj*°(u,vit)  g^^*°\t,o)dt 


(37)-  ♦ K?‘^(u,v;x)g^°»^^(— ,T)dT  + / Jcj (u , v ;t  )g  ^ ^ ^ — ,t  )dx 

o /7  o " /7 

h//2  h//2  , . , . 

♦ / J K^’^(u,v;t,x)g'»^(t,x)dtdx  . 

X=0  tsx 


By  (6) 


K\’°(u,v;t)  = ?[,((u-t)\(»(— ,t  ,u)) 


and  with 


/ . 1-1  for  o < u < t < h//T  I 

I ♦(ib.t.u)  H . ■ ■ ■ 


o foro<t<u<h//7 


K^^’°(u,v;t)  = JSjCCt-u)^)  = (t-u)^  - (t-v)^ 

K^’^(u,v;x)  = ^(^((v-x)^)  Ho 

kV^(u,v;x)  = ?(,((u-— )(v-x)^°)  H o 

(u ,v  ;t  ,x ) = ( (u-t  )'4i (— ,t  ,u)  • (v-x  )°4>(o,x  ,v) ) 

n . i ^ 

- ?{3((t-u)^(v-x)^°)  = ((t-u)^  - (t-v)_j^)  (v-x  ° . 

Similar  arguments  like  in  the  proof  of  Theorem  1 show  that  (here 
:=  {(u,v)l  o < V < u<^>,  :=  T^«l 


w'*  Oi-v 

f Kj«°(.,-;t)g^2’°^t,o)dtll^^ 
o * h 


i h^ll  / K5’°(.,-,t)dtIlL^(^^:,j  . llg^^’°^•,o)llL^^o^h_J 


(39) 


h//2 


< h^'^in'^(K5>°(.,.;t))"dtIlL  (T  ==)  • lo,^5-l 

" — I o 2 1 ^ yo 


/2 


{40)  i 


II//  Kj’^(-,-it,T)6^^’^^t,T)dtdTl| 


L_(V=) 


< h^ll//  kJ’^C- ,•  ;t,T)dtdT|l  • Ilg^^’^^lL  (T  :=) 

- «»n 


and  so  on . Now 


K5’°(u,v;t)  = ^ 


^ _ (u-v)  < o for  o<v$u<t<  1//2 
-Ct-v)<o  for  o<v<t<u$l/i'^ 
o for  o < t < V < u 


and  so 


1//? 


/’  ’°(u,v,t)dt  = - (u-v)(^  - 

* A /7 


/f 


1//7 


/ (Kj’°(u,v,t))  dt  = (u-v)  (.—  - ^^^) 


n 


That  implies 


(41)  II 


o « 1 


/ K5*°(.,.it)dt||j^  «)  = ^ ’ 

O «0  1 

2 


(42)  II 


1//?  , ^ 

/ (K5’°(-.-Vt)>  dt|K  ^ :sj 

o 2 1 


3/1377 
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L„[-h//2.h//2] 


-h//2 


1//2 


< 11  / k5’°(-  .•  ;tMtll^ 

o “>  1 


(45)-^ 


c{  11e^^’°\-  ,o)ll^  ,llc  ’ »°^llL^[-h//2,o 


]) 


1//2 


= II  / kJ’°(- ,•  ;t)dt||^  • |lg”’°^*  ,o)l|j^  I-h//2,h//2] 

o “>  1 “ 


and 


(46) 


IIJ/  K?’^- it,T)g^^’^\t,T)dtdTll  . 

no  n 


< h^ll  ,*  it.T))  cltdi'll^^^^^^..) 


■ "l„(T  ==  T.  " ) 


T,-.= 


2'"1 


"2'  h h,-' 


Finally  we  have  to  use  the  inverse  transformation  to  (35). 
By  (35) 


a 

7TI 


1 (3_  + i_)  ^ 

^ ^ax  ay'*  av 


/7 


(—  - — ) v 
'ay  ax'  * 


o iff  X 


y. 


and  so,  see  for  instance  (39),  (40) 


r 


1 (|_  + 1^)  f(x,h-x)  , 


(M7)  A 


(2,1),  V 1 ,3  3 , ,3  3 ^ . 


Fitting  together  (38),  (41),  (43),  (45),  (47)  we  find  (33),  whereas 
(39),  (42),  (44),  (46),  (47)  gives  (34).  □ 


If  one  imposes  additional  symmetry  conditions  on  f,  one  qan  construct 
interpolating  functions  satisfying  (28)  and  reproducing  bivariate 
polynomials . 
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